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SECTION 1 











not possible of accomplishment due meinly to lack of time, are 


hereinafter presented. 
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Tannerature, Fuel Flow, and Propeller Biees Angle. In addition, 
meximum engine efficiency is attained at neerly maximum allowable 
turbine speed and temperature, thereby requiring operation in a 
regime having narrow control limitations on certain of the veriables. 

The ultimate requirement for a turbo-prop control system 
is that the pilot be able to select a desired value of ynower, without 
exceeding ay of the design limitations of the power plant components. 

A great amount of research is «ct nresent being conducted in an attenot 
to obtain such a control syste. 

In References 1, 2 and 3,a general discussion of the control 
requirements for turbo-vrop power plants is given. In Ref. 4, the 
edaptability of the electronic differential analyzer to the analysis 
of the requirenents of the control systen is discussed. In Ref. 5, the 
utility of the electronic differential analyzer for obtaining solutions 
to megineering problems is considered, and a complete description of the 
type analyzer used in this investigation is given. 

It was the purmose of this investigetion, in addition to giving 
the authors practice in the use of the electronic analyzer, to carry out 
& Study of the basic requirerents of the control system for a turbo-prop 
power plant. Since the purpose did not include the design of such a 
system, a fundemental anprmach has been attenpted throughout. The method 
followed was to start with an initial, simplified equation representing 
the basic charscteristics of the turbo-prop power plant, and to add to the 
initial system, as time allowed, certain functions which would improve 
the r@sponse characteristics of the coutrol systen. 

It is regretted that due to time limitations, a more complete 


luvestigation could not be made. 
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SECTION III A 














Instantaneous eerodynanic torque absorbed by 


the proneller 


Instantaneous net engine torque (equals 2-& | 
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Development of engine equation: 


Note: Constent density ratio, temperature retio, 


and flicat velocity assumed here and also in 
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the development of the propeller equetion. 
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where, Qa = Qa (te, 4) = Qa (Ae rdinc, B% rhs) 
Q. ~ 0@./dn«), dine + (d Q./ 95), . 
) 


il 


Qa, 

Qa = Q.. + Cp. Ae = oo Cres Ad 2 

Q, fe = Ji + Cp,, dre + Cp, AZ, 
Therefore, 


Tp/it dA AeAeVdt = Qpe - /2- Spa de ~ Gs bs. 


But, 
* ae ¥ 
Ope = Qee TY , aho Qa fx = Me ‘ 


Therefore the propeller equation is: 


Lpy,3 A (rete )/se =¢-Cp,, dae —Cos dp. (2) 
DC WiGaaane)/de-= Cay, Abe * Coy, Jt 
Tp fia d(aewe je = $ ~Cpae 422-6. (2) 
Eliminating q and substituting differences for differentials: 
(Ig + Ip /.2) d(Ac-ne ye + Crs (3-3") 
+(Cpae + Cen, )iAe-Me) = Coy, (ee). (3) 


Equation (3) is the basic equation used in this investigation. 








SECTION III B 
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DATA ON TH: T. G. 100, TURBO-FROP ENGIE 


For purposes of this investigation, operating data on the 
General Electric T. G. 100 turbo-prop engine, as presented in Ref. 6, 
was used exclusively since this was the only data available to the 
authors in an unrestricted form. The fact that the T. G. 100 is 
essentially a constant speed engine, does not reduce its value as 
representative of a turbo-prop engine, but does simplify the control 
probl an to some extent. Ina variable speed engine, a new engine speed 
accompanies each new power plant output called for, in order to operate 
at maximum system efficiency. This affects the control system problem 
in that this new engine speed would have to be called for,or scheduled. 

For this investigation, a flight velocity of 300 MPH was 
assumed for Step 1; end for Step II, veloetties of 100, 300, and 500 
MPH were assumed; all at a fuel flow rate of 1400 lbs/hr and at sea 
level. In order to evaluate the coefficients for the equation repre- 
senting the power plant, the values of total shaft horsepower listed 
in Table I were obtained from Figure OA of Section 24, Ref. 6, repro- 
duced in this report as Figure 2. The friction loss was obtained from 


Figure 25 of Section 23, Ref. 6, and subtracted from the total shaft 








horsepower thus obtaining the “ngine Shaft horsepover (SM). This 


im 


was tiuen converted to engine torque ( 


335,000 
i = ee a 
2177 (Ne , RPM) - A plot of Engine 
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Ae ft.tbs=SHe 


( Xe ) versus Engine Speed Ca) with fuel flow rate ( w a. as par 


n made for 100, 300 and 500 IiFH flight velocities in Figures 
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SHCTION III C 











engine power, and engine RFt! vere selected at sea level, and and 

Co computed. Entering Figure 6 with these valuss of ¥ and Cp » points 
indicating 3 and propeller efficiency were determined. Plotting a 
nusiber of these points for speeds of 100, 300 and S00 MH indicated that 
at 300 LPH the points were to the left of maximum efficiency, while for 
500 PH, the points were to the right. This phenomena was caused by too 
large a variation in + due to the fact that an initial basic premise 
was constant wie eo riona velocity for steady-state operation. 

This loss of propeller efficiency was disregarded because, for the pur- 
pose of this analysis, only adequate power absorption of the propeller 
is important, not power output. Since this situation does exist, however, 
it is implied that evon though an engine of itself gives maximum econony 
at full rated RPM, a simaltmeous schedule of fuel flow and RFN instead 
of fuel flow only, is necessary to comvensate for the loss in efficiency 
Of the propeller. 

Having decided that the propeller was acceptable, the arguments 
for Table II were obtained from Figure 6 by entering this Figure with 
end reading off Co for verious 2 » Fron this table, Figures 7 through 
Figure 12 were constructed for the three speeds selected for analysis. 
These curves are plots of <a vs B With Me as parameter and oe VS Sr_E 
wi th 8 as parameter. The purpose of these curves is simply to enable 
determination of Cog and Cong for each operating condition analyzed 60 
that these coefficients might be inserted in the general differential 
equation representing the system. Ce. end Cone oF Slopes at the par- 
ticuler operating condition, being equal to 0a/:/98), 208 (4 //d.n6), , 
respectively. The coefficients for flight speeds of 100, 300, and 500 


MPH at sea level and We of 1400 lbs/hr are tabulated in Table III. 
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One other quantity needed is the moment of inertia of the 


propeller. whether of the solid aluminum alloy or hollow steel type, 





the result was essentially the sme. One method of approach is illus - 


trated in Ref. 8, and is applicable to aluminum alloy blades. The 


estimated Ip the method of Ref. 8 was 70.2 slug feet”. 
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SECTION III & 





The basic equation representing tne power plant is equation 
(53) (see Section III A above). P 
s ty 7 ¥ = _ 
Fg + Ss] acta) + CpaGAse") + ene Cene ] te ee) 


= Ce, (we-5") , (3) 


For purposes of analysis, a new veriable is introduced, i.e.; 
y 9 
@ = J (Agr te )dt, 


@ 


* 
@= (Ne =e), 


o@ 


S) 
Substituting 6 in equation (3), the basic equation becomes: 


[ter Teg] e + Cp) t[Cpae-Cen,] 6 
= Cows (Ws -), (4) 


If it is now assumed that the propeller pitch control chenges 


(Le ~Ae) 


F 


vropeller pitch at a rate vrovortional to the engine offspeed, tha 


AA/e = k, (Ae Ae), | 
we = 3 + ky { Ae Re) dt = 3" + k; @) | 

t : 
AiB* = hf (end dt= ke (5) 
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J ox 


where Ky is the governor proportionality constant. This type of 
control on proveller pitch will be termed “integral control on @ :. 
Substituting this value of (3 -A* ) in equation (4), the 


folloving second order differential equation is obtained: 
[te + ls] & + [age Cag] & +UCps 4] 2 
x 
= Ch (ts ide (6) 

Analysis of the engine and propeller curves, Figure 3 through 
Figure 12, indicates that the coefficients af equation (6) remain prac- 
tically constant over a reasonable range of operation, the inertia term 
netureally remaining constent,. 

If the coefficients are considered as constant, equation (6) 
is conparable to the standard single degree of freedom vibratory equation 


a,e + 2,6+4,96 = F(t), 


where A, = [Te i Tp/:a] is the inertia or mass term, 
uw, =[Cpr- 7 Cons) is the damping tern, 


and Ao = [ Cp. k,] is the spring constant term. 


The known procedure for the ealysis of this type equation can now be 


followed. 


(A) Step l. 
As indicated in Ref. 3, the natural camping term (Cpag ~ Cay, ) 
of the systemis insufficient for adequate response. 


The natural damping was therefore augmented by requiring the 


prop pitch control to change pitch at a rate proportional to the 


‘15 








derivative of engine offspeed, in addition to the engine offsneed. 


_ k, (te De) + ky Ve =e) ) 


A B/ Jz 
BA 3" 


c a *€ 
ky [ tetedelt ik (Ache) . 


k,@ +°k, 6. (Z) 


The basic equation after substitution then becomes 


[Ze + Fof]@ +[Ceae ~ Con, +p] + Gh 
= Cam. (5 _ Ws) ° (g) 


for purposes of analysis, it is convenient to consider the 


ii 


BB 


right side of equation (8) as a step input of engine torque (a. ) 
resulting from a step input of fuel flow rate (Ws; - al ); end then to 
determine the response charecteristics of the system due to such a step 
input. It is considered that this is a conservetive approach, since 
the effect of an actual change in fuel flow rate, involving necessarily 
a finite rate of change, would represent a less stringent condition. 


The basic equation, equation (8), can now be expressed as 


(Te + Tp/,.] © *[Cene~Cere +Cpghul@ + Coa k, @ 


= Gur 3 


Les + Lp/ia] is the inertia term, 


Nt} 


whe re ay 
a, = [ Gpnk . Con, + Co, A, | is the augmented damping term, 


Ap = [ Crs k, ] is the spring constant tern. 
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Various values of damping retio { - = a, /2 Jaa) were next 
assumed; und for cach value of damming ratio, various values of the 
natural undemmed frequency (W,. =(20/a,) of the system were assumed. 
Since the value of the inertia term (a, ) is knowm end remains constant, 
the value of the damping term (s, ) and the spring constant tem (a, ) 
could be determined for each combination of 4 ed Wy chosen: i.e. 


ae fe toe = Ce 4.55), 


aA 
bos 5S w 
= Ay Wr a a i ) 


a, = a§ /a,a. . 


For this part of the investigation, a cruising operating 
condition at sea level was assumed, with a flight velocity ia of 300 
MPH od fuel flav rete (Ww) of 1400 lbs/hr. For this condition, the 
values of Png , re , and “pa had previously been determined and the 
values listed in Table IiI. It was then possible to evaluate the pro- 


portionality constants ky and k,, since 


koa o/ Cp g 


_ Ai -C rE Cerne 
ky 
£8 : 
For completeness of information, the velue of the natural 
demred frequency (MW) of the systen for each combi: tion of 4 and @,, was 


determined by W = W, 1- 42 ° 


The values of a,, 8, Kk), ad k, for the combinations of { ond 


G), chosen are given in Table IV. The reciprocals of the coefficients 


are also given since they ere used in the enalyzer set up. 








All coefficients having been determined, the differential 
analyzer was set up. The diagram of the analyzer circuit is given in 
Fisure 18. Amplifiers A, thro ugh A, are used to solve the basic 
equation (equation (9); while emplifiers hy and AR were added to solve 
the engine equation (equation (1)) for q, in order to record this variable. 
The step input of torque due to change in fuel flow rate (Fu, was applied 
as a battery voltage across the inputs to anplifiers hy and y when the 
switch 5) was closed. An example of the record of a typical step input 
(Fug ) on the Brush recorder is given in Figure 19. 

In order to stay within the voltage limitations of the ampli- 

fiers (approximately 0.5 volts to 100 volts), it was necessary to change 
certain input and feedback amplifier resistances by factors of ten. 
This resulted in some of the variables recorded being different from the 
actual values by a multiplication factor, the factor being indicated in 
Ficure 18. .. correction for this factor was made when the recorded data 
was — 


x Q * 
Tae outputs of amplifiers A,, (te-4e ); A, (A); A, (B78 ); 


a 


and Ag» (q) were recorded. 
For each combination of % andw&,, the proper coefficients, as 
given in Table IV, were inserted in the analyzer in the proper position, 
as indicated in Figure 18. 
A step input, $u,, was then applied for each combination and 
Brush oscillograph traces taken of the response of the variables (te -1s), 
(q), (2 -2"), and (@). The recorded oscillograph traces are shown in 


Figure 14 through Figure 17, respectively. 
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(B) Step 2. 

AS will be seen from the results of this step, the response 
of the control system varies over a wide range with changes in flight 
velocity ifthe type of propeller pjtch control described in Step 1 is 
enployed. If for instance, the proportionality constants ky and ky 
are chosen to give proper response at 300 FH, then these same constants 
coupled with the coefficients corresponding to a different speed cause 
less damping et lower flight velocities and more damping at higher 
velocities. This change in response sith change in speed is illustrated 
by the 0% D.C. curve of Figure 20. (% D.C. refers to the percent deriv- 
ative control used end is defined later on in this section). It would 
be desirsble if the control system could provide approximately constant 
optimum response over the complete range of flight velocity. To 
approach this result we will therefore require the propeller pitch 
control to change pitch at a rate proportional to the second derivative 


of the engine offspeed in addition to the requirements of Step 1. Thus 
+ e * @¢ 
A 4/44 k, (Are she) + ky (tre ~Le) + Ky (Le ~/e) ) 


i) 


Mi 


, t ei ky (a ag) + ky (ene) 
B-8 ki S@rewte)dt + ka e~e 3 ) 


B-B = k, e + ka @ + K3@, 


where K., is enother propeller pitch control proportionality constent. 


This added term can be called “anticipatory"™ or "derivative control on 


B. 
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Again substituting in equation (3), the basic equation for Step 2 


pecomes 


[Ze + tp/:4 + Cog kJ Q + | CG ae Cw k, Je 
ap CP oO f ws ) 
wets ay= [Te tte 4+ Ces k;] is the inertia tern, 
a, = [Crn, + Come + Cp, k| is the damping term, 


ne [ Cp, an | is the spring constent tern. 


It is seen that the addition of derivative control ong cen 
be considered bs en addition to the inertia term (a,) of equation (10), 
thus giving a new “effective” inertia. The basic equation remains in 
the form of the stamdard vibratory equation and can be analyzed in a 
manner similar to that used in Step l. 

An undamped natural frequency (o ) of 5 rad/sec was now arbi- 
trarily chosen as representative of this type power plant at an operating 
condition of 300 MPH, end a damping ratio (4) of 0.6 was chosen as an 
Optimum overall value at this speed. Since the derivative control on 
4 cen be considered to change the inertia term (a,) of equation (10), 
the rrocedure folloved was to increase the value of the inertia term 
by convenient percentages, i.e. 100%, 200%, and 300%. This percent 
increase was termed percent derivative control, amd is abbreviated here- 
after as "% D.U." 

in order to cover the complete ranse of flight velocities, 
three representative operating conditious were chosen, i.e. low velocity 


of 100 MPH, cruising velocity of 30 MH, and high velocity of 500 XFH; 


a0 








all st sea level end a fuel flow rate of 1400 lbs/hr. ‘The determina- 


tion of the velues of the coefficients Ces ; Cpa, » and wy for these 
operatiig caditions has already been described, and the values ar 


siven in Table III. 


By naintaining 











Cp ke 


? C pre = Ore + Cpe ka 


- [it Roc] [Te +2] 














From the oscillograph of the four variables recorded, the 


maximum and steady state deviations in millimeters, slope in the case 
of B ; and frecuermy and time to damp were obtained. These measure- 
ments in millimeters were translated into voltages by applying the 
calibration data and attenuator settings on the recorders, and also 

the multiplication factors made necessary by the voltage limitations 
of the analyzer amplifiers. These voltage magnitudes were tren eqal 
quantitatively to the value of the variable measured. The values of 
frequency and time to demp were difficult to obtain from the traces 
end therefore the absolute values given can be considered only approxi-~ 
mate. Values of frequency not listed in the tables were impossible to 
messure from the traces. «an examle of the method used in the reduction 


of data is given in Section VII, Sample Calculations, 


Step 1. Results 
Step 1 was based on an original steady stete torque of about 
65% power, The step input of torque was 63.5 ft. lbs. or equivalent to 
roughly a 10% change in power based on the original steady state power. 


This corresponds to about a 6.5% change based on maximum power. Since 
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the coefficients of the differential equation for the system change 
only slightly with change in output of the engine, the equation can 

be assumed linear for all values of torque at one airspeed without 

too such error. For instance, if a jump from 65% to 100% torque were 
desirable, the results tabulated could be multiplied by 100-65/6.5=5.4 
to give fairly good approximations of response. 

In this step, it was assumed that @ was a function of the 
integral of the offspeed ed the offspeed. The traces taken are pre- 
sented in Figure 14 through Figure 17, and are arranged to simw the 
effect on response of changes in 9 end wo: Only relative amd not 
quantitative comparisons of maximum and steady state values of Oe am 
@ shorm on the traces may be made, due to the necessity of having dif- 
ferent attenuation end calibration values set on the recorders. This 
was necessary in order to get traces of reasonable size. The reduced 
data is given in Table IV, and Figure 25 was drawn from the tabulated 
values. 

Figure 25 shows the effect on maximum response values of q, 
ot: f ana 4 of varying damping ratio (4) and natural undamped fre- 
quency (W). It is seen that increasing W, results in lower maximum 
values of tarque and engine offsveed, end higher maximum values or 4 
end af - Increasing ¥ , in general, decreases the maximum values of 


dt 


the variables plotted except for AZ in which case the reverse is true. 


At 
In an actial control system, there would naturally be limits 
AZ ng 4 
on the maximum allowable or obtainable re and ce: These limitations 


would partially define the design problem which could be solved by 


selecting proportionality constents (k, and k,) to obtain the necessary 


‘4 ndW, for adequate response. 
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d22 
dt+ 
due to the inability of measuring this quantity in all cases from the 


wes not caasidered previously in this investigation 





AB traces. Several apnvroximate values were obteined from the % = 50 


(42 [sect 


At Waid 


traces; 








in speed, due to maintaining constant fuel flow instead of constant 
torque throughout the sneed range. The comments mede for Step 1 


recarding lineer assumptions apply here also. 


In this step, derivative control on 4 was added to the p 


vious requirements of the propeller pitcn control to col 











of we could not be measured. However, a plot of the change in 8 at 
a 


this large initial $4 is siven in Figire 27. By inference, it may 

















be theorized from Figure 27 that the lower the speed and the greater the 
amount of derivative control, the more extreme the demands or require- 
ments on the propeller cmtrol mechenism will be. This narticuler 


characteristic of response using derivative control sould definitely 


















be investigated further. Cne way measurements of this phenomena could 
be easily accomplished with the cifferential malyzer is by changing 


the time base. 


It is realized that the turbo prop control problem is very 
\ 
camplex end has not been solved to anyone's comclete satisfaction to 
date, as is evidenced by tre fact that the first production turbo prop 
systen in the United States has only recantly been installed in the 
Navy's XPSY. This investigation has been based on possibly rudimentery 


considerations but a sincere attempt has been made to be as general as 


possible. I¢ is realized that the investigation is not complete end | 





that certain other approaches are desireble, even though not possible 
‘to investigate, due to lack of time. A discussion of some of these 
other consideratias follows. 

A comlete system analysis would necesssrily consider the 
Pedi ining variable not touched upon in this investigation, nee 


i 
7 


turbine bucket termerature. Such a system, as visualized by the authors, 





wolld introdwce this variable as a limit to fuel flow rete. d‘hat the 
functimal relaticnship betveen fuel flow rate and temperature should 


be or might be has not been investigated. Hovever, it is apparent that 


rete of chenge of fuel flow rate must be so limited, regerdless of 


99 

















throttle quadrant movement, to values such that temerature limita- 
tions will not be exceeded. In any event, the result would be less 
stringent than a step input of fuel flow, as is being used in this 
investigation, and results nerein obtained Should therefore be con- 
servative. 


One basic assumption was to consider rotational speed 



















governed exclusively by propeller blade angle. One problem envisioned 
in such an assumption would be in the case of simultaneous scheduling _ 
of RPM and fuel flow, which is in reality necessary as mentioned in 

the section on the derivetion of oropeller data (III C). In this case, 


an increase in setting of fuel flow and RPM simuilteneously would cause 


am initial decrease in blade mele due to the underspee® comition, 





resulting in an initially lower rather than higher torque amd probably | 
also in mare extreme oscillations. A discussion of this is civen in 


Ra “@ 2 
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Propeller blede mgle control of RPM is not, however, the 


only approach possible, although seemingly the most feasible and practi- 





cable, «An alternate approach as mentioned in Ref. (2), would cause fuel 





) 


flow ¢ to bo a function of RPM only, with propeller blade mele ca ntro 1] 


i 1 5 | Y 1 rn 
i Ling 
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the remaining variables, 





No time lags, lost motion, or finite limitations on any of 





variables were considered. This would be of importance in a physical 


System where, for instance, a dé were required instentanesously 








control mechanism or in the reduction gearing. and also the probability 


) 
that the fuel pump response would be non-linear and asymptotic to both 
Origimil end new steady stcte velues of fuel flow. 


Such phenomena would heve the effect of non-linearizi 


system with attendant effect on response. It is con 
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ng¢ term would then 


E +] §) +0] 6] ; 


It is not denied that the dssign of a workeble systen 
utilizing these principles would de difficult. Differentiation, inte- 


eration, addition, subtraction, multiplication, and division of variables 











ae Lan a 


woild have to be accomplished in the control instellation. Yrresentation 
of tiis informtior, therefore, is intended to show nerely whet might be 


accoiplished theoretically, uSing differential analyzer methods. 

No effort Las been made to indicate any physical hardwere 
a 
ai 


p 
ops 


ich would be necessary to set up an actual control System on ony 


tue premises investigated or discussed, since this investigation ¢ 








fe 


ith investigation, not with design. It is considered, however, thet 


™ ctronic differential analyzer non-linear ca.ponents could be us 


would be operated on, resulting in the formulation of correction signe 









S visualized by the authors, one such system would operate with 





control an 








blade angle or fuel flow, or voth,as deterrined to be desirable. 


Since varieble demping gives the best response, tuis might be incor-= 


porated in the propeller control. A movement of the throttle co 











SBCTION V 
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derivative, a consistent systan response at all operzting 
conditions may be approached by the addition of the vrover 


amount cf derivetive control. 


bh, 
i, 
re 


The propeller pitch rate end acceler tin should b 


orecticel veluos to 
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